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We present measurements of the top quark mass using mT2, a variable related to the transverse 
mass in events with two missing particles. We use the template method applied to tt dilepton 
events produced in pp collisions at Fermilab's Tovatron Collider and collected by the CDF detector. 
Prom a data sample corresponding to an integrated luminosity of 3.4 fb"'^, we select 236 tt candidate 
events. Using the mT2 distribution, wc measure the top quark mass to be Mtop = 168.0 j]* * (stat.)± 
2.9 (syst.) GeV/c^. By combining mT2 with the reconstructed top quark mass distributions based 
on a neutrino weighting method, we measure Mtop = 169.3 ± 2.7 (stat.) ± 3.2 (syst.) GeV/c^ This 
is the first application of the mT2 variable in a mass measurement at a hadron collider. 
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I. INTRODUCTION 

Models in numerous, well-motivated theoretical frame- 
works make predictions for new phenomena at hadron 
colliders such as the Tevatron and the Large Hadron Col- 
lider (LHC) [l|,[2. Within each framework, one can con- 
struct a number of qualitatively different models consis- 
tent with data. Thus, when discoveries are made at a 
hadron collider, we face the inverse problem of how one 
maps back to the underlying theory responsible for the 
new phenomena [U, Q • A potentially powerful observable 
to discriminate among models and to extract the mass of 
new particles, when the new phenomenon produces a pair 
of new particles with large missing energy signatures, is 
the mx2 variable [3, Q ■ The r7ix2 variable is based on 
transverse mass in events with two missing particles. 

The top quark is the heaviest known elementary par- 
ticle with a mass approximately 40 times larger than the 
mass of its isospin partner, the bottom quark (6). The 
large top quark mass (Mtop) produces significant contri- 
butions to electroweak radiative corrections. Therefore, 
top quark mass measurements are important tests of the 
standard model and provide constraints on the Higgs bo- 
son mass. In the dilepton channel, tt pair production 
is followed by the decay of each top quark to a W bo- 
son and a b quark where both W bosons then decay to 
charged leptons (e or fi) and neutrinos. Events in this 
channel thus contain two leptons, two b quark jets, and 
two undetected neutrinos. The measurement of Mtop us- 
ing complementary techniques tests and improves our un- 
derstanding of this important parameter in the standard 
model 

In this letter, we present the first measurement of the 
mass of the top quark using the 771x2 distribution with tt 
events in the dilepton channel [7] . We use this channel be- 
cause it has decay products similar to possible new phe- 
nomena where undetected particles are created. We com- 
pare this method with two others that were previously 
used: the reconstructed top quark mass using the neu- 
trino weighting algorithm (777.^^^) 8, ^ and the scalar 
sum of transverse energies of jets, leptons, and missing 
transverse energy {^t) LlQ] in the event (Ht) lUl]. We 



also measure the top quark mass using pairs of observ- 
ables ((7^x2, mf^^) and {mf^-^, Ht)) simultaneously. 



II. THE mT2 VARIABLE 

Many models contain heavy, strongly- interacting parti- 
cles with the same conserved charge or parity that result 
in weakly-interacting, stable particles in the final state. 
A hadron collider would pair-produce these colored par- 
ticles, which then decay into standard model particles 
along with a pair of undetectable weakly interacting par- 
ticles, so that the generic experimental signature is large 
missing transverse momentum accompanied by multiple 
energetic jets and leptons • In this final state, we can 
define 7Mx2 as 



771x2 [m^. 



mm 

(1) (2) 

Pt 'Pt 



max[7;7.T("^i™is; Pt'')) »77T("imOTs; PtO] 



(2)n 



(1) 



where rriT, the transverse mass of each parent particle, 
is defined as 



rriTim 



. mvis \ 
invis 5 Pt J 



2{E^'''E'^-"'' - p^"* • ptp""'). (2) 



Here invis and vis represent the individual unde- 
tected (invisible) and detected (visible) particles, respec- 
tively, p^'* and p^^ are transverse momenta of two in- 
visible particles and m.invis is the mass of the invisible 
particle. The minimization is performed with the con- 



straint p^^^ -I- p^^ 



Pjn , where the magnitude 01 



prp " IS constrained to the missing transverse momen- 
tum. 

The quantity 777x2 represents a lower bound on the 
mass of the parent particle. Using the 777x2 distribution, 
we can extract the mass of this parent particle [l^ in a 
similar way to the precise measurement of the W boson 
mass (isj where an event contains one charged lepton (e 
or jj,) and a neutrino, with the latter not being detected. 
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III. EXPERIMENT AND DATA 

We use a sample of tt candidates in the dilepton chan- 
nel, corresponding to 3.4 fb~^ of proton-antiproton colli- 
sions at -^/i = 1.96 TeV, collected using the CDF II detec- 
tor [3]. This is a general-purpose detector designed to 
study pp collisions at the Fermilab Tevatron. A charged- 
particle tracking system, consisting of a silicon microstrip 
tracker and a drift chamber, is immersed in a 1.4 T mag- 
netic field. Electromagnetic and hadronic calorimeters 
surround the tracking system and measure particle en- 
ergies. Drift chambers and scintillators, located outside 
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TABLE I: Expected and observed numbers of signal and 
background events assuming tt production cross section crjj 
— 6.7±0.8 pb and Mtop ~ 175 GeV/c^. Uncertainties quoted 
capture the uncertainties on the theoretical cross section, the 
statistics of data in the Z mass window, the jet energy scale, 
the luminosity, the fake rates, and the statistics of the MC 
samples. 
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FIG. 1: The mT2 distributions from tt dilepton Monte Carlo 
events that pass the selection criteria for three input values of 
the top quark mass. Each distribution is normalized to have 
unit area. 



the calorimeters, detect muon candidates. 

We select events consistent with the tt dilepton decay 
topology. We require two oppositely charged lepton can- 
didates with pt > 20 GeV/c with one isolated [lal lepton 
candidate in the central region (|r/| < 1) of the detector, 
and another isolated or non-isolated lepton candidate in 
the central region, or isolated electron candidate in the 
forward region (1.0 < \ri\ < 2.0). We also require ex- 
ceeding 25 GeV, and at least two jets with Et > 15 GeV 
and |?7| < 2.5 [Toj- To further reject backgrounds, we 
request Ht > 200 GeV. We also require the variables 
of interest to be consistent with the top quark hypothe- 
sis by demanding 20 GeV/c^ < mT2 < 300 GeV/c^ and 
100 GeV/c^ < m^^^ < 350 GeV/c^. The criteria select 
236 tt candidate events. 

The primary sources of background production are 
Drell-Yan, diboson, and QCD multijet events. We es- 
timate the rate of the Drell-Yan events with a calcula- 
tion based on simulated events using the ALPGEN fl^ 
v2.10 Monte Carlo (MC) generator and the rate of dibo- 
son events with a PYTHIA [13] v6.216 calculation. For the 
Drell-Yan Z-|-jets process, we normalize the MC sample 
by matching the number of Z events predicted and ob- 
served in the Z mass region between 76 GeV/c^ and 106 
GeV/c^ . We use data to estimate the rate of background 
events from QCD multijet production where an event has 
one real lepton and one of the jets misidentified as an- 
other lepton (fake). In measuring the top quark mass, 
we divide the tt candidate sample into events with and 
without secondary vertex b tags [18], which have very dif- 
ferent purity. We only attempt to b tags the two highest 
Et jets. Table U summarizes the composition of back- 
ground events and the expected numbers of tt and back- 
ground events. We estimate the tt signal event rates us- 
ing PYTHIA v6.216 with CTEQ5l parton distribution 
functions at leading order with a full detector simula- 
tion d^. 

To calculate tot2 of a tt dilepton event we first 
identify all possible configurations corresponding to dif- 



ferent assignments of jets to b quarks and combinations 
of quarks and leptons. The two most energetic jets in an 
event are considered to have originated from the b quarks. 
For each configuration, we calculate the transverse mass 
of each top quark (t — > bli^) using Eq. [2j 

ruT = ,Jml +ml + 2{E!^E!^ - • p^), (3) 

where mt,; and p^ denote the invariant mass and trans- 
verse momentum of the bottom-quark jets and charged 
lepton (bl) system, and are the mass and trans- 
verse momentum of the neutrino, and E"^' and Elf are 
the transverse energies of the bl system and neutrino: 

Et = ^J\Pt\^+K and Ef = ^\plf\^+ml. (4) 

We then calculate r7iT2 using Eq. [1] with the assump- 
tion m,y = 0, and for all possible parton assignments. 
We select the smallest value for each event. Figure [T] 
shows simulated mT2 distributions for various top quark 
masses for the combined non 6-tagged and 6-tagged sam- 
ple, which demonstrates that tot2 is sensitive to Mtop, 
and thus can be used to measure it. 



IV. MASS FIT 

We estimate the probability density functions (PDFs) 
of signals and background using the kernel density es- 
timation (KDE) ^21, 22] that constructs the PDF with- 
out any assumption of a functional form. For the mass 
measurement with two observables, we use the two di- 
mensional KDE that accounts for the correlation be- 
tween the two observables. First, at discrete values of 
Mtop from 130 GeV/c^ to 220 GeV/c^ with increments 
of 0.5 GeV/c^ in the region immediately above and below 
175 GeV/c^ to 5 GeV/c^ near the extreme mass values, 
we estimate the PDFs for the observables from 76 ti MC 
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samples. Each sample consists of 0.6M to 4.8M gener- 
ated events, with IM events corresponding to a luminos- 
ity of 150 fb^^, assuming a tt cross section of 6.7 pb [23l |. 
We smooth and interpolate the MC distributions to find 
PDFs for arbitrary values of Mtop using the local poly- 
nomial smoothing method [24]. We fit the distributions 
of the observables in the data to the signal and back- 
ground PDFs in an unbinned extended maximum likeli- 
hood fit [25} where we minimize the negative logarithm of 
the likelihood using minuit [11]. The likelihood is built 
for the 6-tagged and non ^-tagged categories separately 
and then combined by multiplying the two categories. 
We find the statistical uncertainty on Mtop by search- 
ing for the points where the negative logarithm of the 
likelihood minimized with respect to all other parame- 
ters deviates by 0.5 units from the minimum. Ref. (22j 
provides detailed information about this technique. 

We test the mass fit procedures using 3,000 pseudo- 
experiments for each of 14 different top quark masses 
ranging from 159 GeV/c^ to 185 GeV/c^ with almost 
2 GeV/c^ step size. In each experiment, we select the 
numbers of background events from a Poisson distribu- 
tion with a mean equal to the expected numbers of back- 
ground events in the sample and the numbers of signal 
events from a Poisson distribution with a mean equal 
to the expected numbers of signal events assuming a tt 
pair production cross section of 6.7 pb. The distribu- 
tions of the average mass residual (deviation from the 
input top mass) and the width of the pull (the ratio of 
the residual to the uncertainty reported by minuit) for 
simulated experiments show that the measured top quark 
mass IS on average 0.26±0.10 GeV/c^ lower than the true 
top quark mass and has no dependence on Mtop in the 
mT2 measurements. We correct the measurement for this 
bias. No such bias is observed with the combined (mT2, 
j^NWA-j jneasurement. In all cases, the fit on average 
correctly estimates the statistical uncertainties, based on 
the pull width distribution being consistent with unity. 
For Mtop =175 GeV/c^, we expect the statistical uncer- 
tainties on Mtop to be 4.0 GeV/c^ with mT2, 3.4 GeV/c^ 
with 5.4 GeV/c^ with Ht, 2.9 GeV/c^ with (mT2, 

mfWA) combined, and 3.2 GeV/c^ with (mf^A^ 
combined. 



V. SYSTEMATIC UNCERTAINTIES 

We examine a variety of systematic effects that could 
affect the measurement by comparing MC simulated ex- 
periments in which we vary relevant parameters within 
their systematic uncertainties. The dominant source 
of systematic uncertainty is the light quark jet energy 
scale (JES) [13]. We vary JES parameters within their 
uncertainties in both signal and background MC gener- 
ated events and interpret the shifts as uncertainties. The 
&-jet energy scale systematic uncertainty arising from our 
modeling of b fragmentation, b hadron branching frac- 
tions, and calorimeter response captures the additional 



uncertainty not taken into account in the light quark jet 
energy scale. The uncertainty arising from the choice of 
MC generator is estimated by comparing MC simulated 
experiments generated with pythia and herwig [28| . 
We estimate the systematic uncertainty due to model- 
ing of initial-state gluon radiation (ISR) and final-state 
gluon radiation (FSR) by extrapolating uncertainties in 
the pt of Drell-Yan events to the ti mass region [29] . We 
estimate the systematic uncertainty due to parton distri- 
bution functions by varying the independent eigenvectors 
of the CTEQ6m tSOl] parton distribution functions, vary- 
ing Aqcd, and comparing CTEQ5l [19] with mrst72 [31| 
parton distribution functions. In estimating the system- 
atic uncertainty associated with uncertainties in the top 
quark production mechanism, we vary the fraction of top 
quarks produced by gluon-gluon annihilation from 6% to 
20%, corresponding to the one standard deviation up- 
per bound on the gluon fusion fraction js^]. We esti- 
mate systematic uncertainties due to the lepton energy 
and momentum scales by propagating shifts in electron 
energy and muon momentum scales within their uncer- 
tainties. Background shape systematic uncertainties ac- 
count for the variation of the background composition. 
In addition, we change the shape of the Drell-Yan back- 
ground sample according to the difference in the missing 
energy distribution observed in data and simulation, and 
the shape of the QCD multijet model. We estimate the 
multiple hadron interaction systematic uncertainties to 
account for the fact that the average number of interac- 
tions in our MC samples are not equal to the number 
observed in the data. We extract the mass dependence 
on the number of interactions in MC pseudo-experiments 
by dividing our MC samples into subsamples with differ- 
ent number of interactions. We then multiply the slope 
of the result by the difference in the number of interac- 
tions between MC events and data and treat that as a 
systematic uncertainty. 

It has been suggested that color reconnection (CR) ef- 
fects could cause a bias in the top quark mass measure- 
ment and interpretations at the level of 0.5 GeV/c^ [33| . 
We estimate uncertainties arising from CR effects using 
the PYTHIA 6.4 MC generator, which includes CR effects 
and other new features in modeling the underlying event, 
initial and final state radiation, and parton showering. 
We generate two MC samples, one using tune A [sj], 
which is very similar to the tune for CDF nominal MC 
generations, the other using Acr [33], which includes CR 
into the tune A. We take the difference in the extracted 
mass between these two MC samples as a systematic un- 
certainty. We measure the difference to be 0.6 GeV for 
(tox2: m^^^) combined, and 0.7 GeV for tox2 alone. As 
a cross check, we generate two other MC samples, one 
using tune So [l^ and the other using NOCR [s^, which 
include all of the new features with and without CR. We 
find a similar mass difference between the two samples. 

Table Hi] summarizes the sources and estimates of 
systematic uncertainties. The total systematic uncer- 
tainties, adding them in quadrature, are 2.9 GeV/c^ 
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TABLE II: Estimated statistical (Mtop =175 GeV/c^), systematic, and total uncertainties in GeV/c^. 





771T2 




Ht 




(m^^WA H ) 


Statistical 


4.U 


Q A 


0.4 


o n 

z.y 


Q O 
O.Z 


Systematic 


Jet energy scale (light quarks) 


Z.O 


o c 
O.O 




O A 


0.4 




Generator 


A O 
U.O 


i.U 


Z.D 


U.O 


i.O 




Parton distribution functions 


U.5 


U.D 


1 o 
i.O 


U.O 


A O 

U.O 




b Jet energy scale 


A O 


A O 


A O 

(J.z 


U.^ 


A O 

U.O 




Background shape 




u.o 


n 7 

U. i 


U.i 


u.o 




Gluon fusion fraction 


n 1 
yj.o 




u.o 


\U. i 


n 1 

U.l 




Initial and final state radiation 


0.6 


0.2 


0.6 


0.3 


0.2 




MC statistics 


0.3 


0.3 


0.5 


0.3 


0.3 




Lepton energy 


0.6 


0.2 


0.7 


0.3 


0.2 




Multiple hadron interaction 


0.2 


0.3 


0.3 


0.3 


0.3 




Color reconnection 


0.7 


0.6 


2.5 


0.6 


0.6 




Total systematic uncertainty 


2.9 


3.8 


5.7 


3.2 


3.8 


Total 


5.0 


5.1 


7.8 


4.3 


5.0 



with mT2, 3.8 GeV/c2 with mf^A^ 5 7 GeV/c^ with 
Ht, 3.2 GeV/c2 with (mT2, mf^A) combined, and 
3.8 GeV/c2 with (mf^A^ jj^) combined. The 
method has a jet energy scale uncertainty significantly 
smaller than mf^^, resulting in the smallest total sys- 
tematic uncertainty. Including both statistical and sys- 
tematic uncertainties, we conclude that 771x2 is one of 
the best observables for the Mtop measurement , compa- 
rable to the measurement using 777,^^^. Using both 777x2 
and mf^-^, we expect to achieve a 10% improvement in 
overall uncertainty over using 777x2 alone. 



VI. RESULTS 

We apply a likelihood fit to the data using observ- 
ables discussed in this article. Figure [5] shows the one- 
dimensional log-likelihoods for 777x2 and (777x2, mf^^) 
combined. Figure [3] shows the distributions of the ob- 
servables used for the Mtop measurements overlaid with 
density estimates using tt signal events with Mtop = 
169 GeV/c^ and the full background model. The fit re- 
sults are summarized in Table IIIII The extracted masses 
are consistent with each other and the statistical uncer- 
tainties are consistent with predictions from MC pseudo- 
experiments. 

In conclusion, we present the top quark mass measure- 
ments in the dilepton channel using 777x2- In 3.4 fb"' of 
CDF data, we measure Mtop using 7?7x2 to be 

Mtop = 168.0 +to (stat.) ± 2.9 (syst.) GeV/c^ 
= 168.0 tl'^o GeV/c^, 



and using both mf and ?77x2 to be 

Mtop = 169.3 ± 2.7 (stat.) ± 3.2 (syst.) GeV/c^ 
= 169.3 ±4.2 GeV/c^. 

This is consistent with the most precise published result 
in this channel from the CDF '35] and DO [36] collabo- 
rations. We expect further improvements in Mtop with 
these variables as CDF accumulates about a factor of 3 
more data during Tevatron Run II. The measurements in 
this article are the first application of the 777x2 variable to 
data, and demonstrate that 7f7-x2 is a powerful observable 
for the mass measurement of the top quark in the dilep- 
ton channel. The methods described in this article will 
be applicable to other measurements at the Tevatron and 
soon at CERN's Large Hadron Collider for discriminat- 
ing new physics models and measuring the mass of heavy 
particles that decay into weakly interacting particles such 
as dark matter candidates. 
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